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ABSTRACT 

New Hi observations of Messier 31 (M31) obtained with the Effelsberg and Green Bank 100-m telescopes 
make it possible to measure the rotation curve of that galaxy out to 35 kpc. Between 20 and 35 kpc, the 
rotation curve is nearly flat at a velocity of ~ 226 km s"' . A model of the mass distribution shows that at the 
last observed velocity point, the minimum dark-to-luminous mass ratio is ~ 0.5 for a total mass of 3.4 x 10^' 
M.Q cit R < 35 kpc. This can be compared to the estimated MW mass of 4.9 X 10" AlgforiJ < 50kpc. 
Subject headings: galaxies: halos — galaxies: fundamental parameter (mass) — galaxies: individual (M3 1) — 
galaxies: kinematics and dynamics galaxies: structure — Local Group 



1. INTRODUCTION 

Over the last 35 years, accurate rotation curves (RC) have 
been derived for a large number of galaxies in the Local Uni- 
verse, using Hi or CO at radio wavelengths, or emission lines 
such as Ha in the optical (see Sofue & Rubin 2001, and ref- 
erences therein). With larger telescopes and new instrumenta- 
tion, RCs have e ven been measured for galaxies at redshifts z 
> 0.5 (Yost et al. 1997^ .7iegler et al..2003^.Flores et al..200l 
Hamfordet al. 2005). 

Surprisingly, the two nearest massive galaxies, the Milky 
Way (MW) and M3 1 (at an adopted distance of 780 kpc from 
McConnachie et al. 2005), have very poorly defined RCs. Our 
position inside the Milky Way's disk makes it very difficult 
to interpret the Hi outside the solar radius, while for M31, 
the gas on the receding side is at the same velocity as MW 
gas and it is difficult to disentangle the two components. As 
importantly, the very large angular size of this galaxy makes it 
a difficult target for synthesis telescopes (like the VLA) which 
are not sensitive to emission from large angular scales. 

In spite of numerous attempts to determine the dynami- 
cal mass of M31, a galaxy rich in globular clusters, satel- 
lites, Hl l regions, planeta ry nebulae and having an extended 
Hi disk (*H artwick & S^gent. 1 974 : Gunn 1975; Rood 1979; 
[van den B erg|ll981t lHodgelll992t Courteau & van den Berg 
[T999t iMerrett et alJ l2003t llbata et al.i i2004<) . mass esti- 
mates are quite far from satisfactory. Indeed, we are 
still uncertain as to whether M31 or the MW is the 
most massive member of the Local Group (lEvans et alJ 
120001) ■ The kinematics of M31 was studied extensively 
about 20 - 30 years ago in Hi (Roberts & Whitehurstlll97£ 
Emersonlf! 976; Cram, Roberts & Whitehurst 1980; Unwir 



1983; Brinks & Shane 1984) and Ha (.Rubin & Ford 197C 
Boulesteix et al. 1987: Kent 1989.) . More recent Hi synthesis 



Laboratoire d'Astrophysique Experimentale (LAE), Observatoire du 
mont Megantic, and Departement de physique, Universite de Montreal, CP. 
6128, Succ. Centre- Ville, Montreal, Qc, Canada H3C 3J7 
Electronic address: Claude. Carignan@Umontreal.ca 

- Observatoire de Paris, section Meudon, GEPI, CRNS-UMR 8111 & Uni- 
versite Paris 7, 5 PI. Janssen, 92195, Meudon, France 

Max-Planck-Institut fur Radioastronomie, Auf dem Hiigel 69, Bonn 
53121, Germany 

* Some of the observations reported here were obtained with the 100 m 
telescope of the Max-Planck-Institut ftir Radioastronomie at Effelsberg 

^ National Radio Astronomy Observatory, P.O. Box 2, Green Bank, WV 
24944, USA. The National Radio Astronomy Observatory is operated by As- 
sociated Universities, Inc., under a cooperative agreement with the National 
Science Foundation. 



work was done bv lBraunllfT99(3) and collaborators using the 
VLA, WSRT and the GBT fThilker et al."2004). 

Modern 21 cm Hi receivers are up to an order-of-magnitude 
more sensitive than those used in the earlier surveys, and with 
100-m telescopes such as those at Effelsberg and Green Bank, 
it is now possible to trace the Hi in M31 to a significantly 
larger radius than before. We have thus begun a program of 
observing 21 cm emission from the outer parts of this galaxy 
at R > 90' along its major axis. To avoid confusion with 
Galactic Hi, only the approaching (SW) side of the galaxy 
was observed. The major-axis of M3 1 has a position angle of 
38° and its central position given in the Nasa/Ipac Extragalac- 
tic Database is ajzooo = 00H2'«44:^3, Snom = +41°16'09". 
Throughout the article, line-of-sight velocities are given in the 
heliocentric rest frame. 

2. OBSERVATIONS 

The first observations were obtained with the Effelsberg 
100-m telescope which has a HPBW = 9'.3. Spectra were 
measured every 4.'5 along the major axis of the S W (approach- 
ing) side, beginning at R = 90', or about 20 kpc from the 
galaxy's center. Spectra were measured modulating between 
the signal band and a comparison band offset 1 .56 MHz, or 
about -300 km s"' . The two-channel HEMT receiver was fol- 
lowed by a 1024 channel autocorrelator split into two bands 
of 512 channels, with a channel spacing of 0.64 km/s. The 
system temperature was < 30 K. Both polarizations were av- 
eraged to reduce the noise. 

The brightest Hi emission in M31 comes from a broad, 
ring-like structure in the inner parts of the galaxy JEmersonI 
[1^76; Unwin 1983; Braun 1990), and it is possible that some 
of this emission might be picked up in a telescope sidelobe 
(^eich et al. 1978) and mask any weak Hi signal from the 
outer parts of the galaxy. Care was thus taken during the Ef- 
felsberg observations to ensure that the main sidelobes in the 
elevation direction were not falling on the galaxy. Most of the 
observations were done at parallactic angles ^ ^5° to -55°, 
which is ^ 90° from the M31 major axis. 

Sidelobes are less of a problem with the 100-m Green Bank 
Telescope (GBT), whose offset optics give it a very clean 
main beam with a first sidelobe level more than 25 dB be- 
low the peak. The GBT's gregorian 21 cm receiver system 
overilluminates the subreflector somewhat, producing a broad 
forward spillover lobe which contains about 4% of the tele- 
scope's total response. However, this sidelobe has a diameter 
> 30° on the sky and its effects on the current observations 
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should be negligible (Lockman & Condon 2005). 

The GBT was used to measure Hi spectra at ten positions 
along the SW major axis of M31 in August and September 
2005. At 21 cm, the GBT has an FWHM angular resolution 
of 9.'1. Spectra were taken with 2.5 kms"' effective velocity 
resolution over a range > 1000 kms"' centered on the M31 
velocity. Frequency switching 'in band' gave good instru- 
mental baselines and good sensitivity. The dual-polarization 
receiver had a system temperature Tsys ^ 18 K. The ten posi- 
tions were observed for times varying from a few minutes to 
more than one hour, with a typical value being 40 minutes. 

Hi emission was detected out to = 45 kpc with the Ef- 
felsberg dish, but here we will discuss only data at < 35 kpc 
where spectra were measured with both Effelsberg and the 
GBT. Moreover, the Hi disk of M3 1 is known to be slightly 
warped (e.g. Newton & Emerson 1977; Henderson 1979). 
Knowing our beamsize and estimates of the warp parameters, 
we have calculated that Hi profiles obtained along a constant 
line of 38° are reliable only up to 35 kpc. 

Within this range the agreement between data from the two 
telescopes is good. For the positions observed with both tele- 
scopes (27.6, 32.7, 34.7 kpc), the mean difference in radial 
velocity is less than 2.5 kms"' , and at these three locations 
the radial velocity is thus given by the average of the two ob- 
servations. The Hi spectra, baseline subtracted, with the best 
fit Gaussian profile superposed, are shown in Figure[2 While 
there is some scaling differences between the profiles coming 
from the two telescopes, the profile shapes are remarkably 
similar (e.g., the two spectra atR = 27.6 kpc). 

3. ROTATION CURVE OF M31 

There are discrepancies between the different determina- 
tions of the RC of M31 in the inner 10 kpc (j^ofue & Kato 
119811) . but there is general agreement for 10 < 7? < 30 kpc. 
One interesting feature of the Hi RC of M3 1 is that it seems 
to be declining more or less regularly from the center all 
the way out to 30 kpc (Braun 1991) . While there are a 
few galaxies with declining RCs over a limited radius range 



TABLE 1 

Hi rotation velocities computed from data IN lUNWIN|jI983D 
FOR THE INNER PARTS (R < 90') AND FROM THE EFFELSBERG AND GBT 
SINGLE DISH OBSERVATIONS FOR THE OUTER PARTS (R > 90')- 
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Notes on columns: (1) and (2) Radius in arcmin and kpc (resp.) for D = 780 
kpc iMcConnachie et al. 2005) (l'= 227 pc) (3): Velocity in km s~ ' , for a 
systemic velocity of —300 km s"' Ide Vaucouleurs et alll99ll) . 

(»Cari2nan & Puche 199^, ICasertano & van Gorkoii? 1199 It 
Honma & Sofue 1997), this is quite unlike what is seen in 
most spirals, and is one of the reasons we wished to obtain 
the velocity information at larger radii to check whether that 
trend continues. 

Tableland Figure |2l give our derived Hi RC. For R < 90' 
(or ~ 20.5 kpc), we have redete rmined the Hi velocity field 
using the early Hi data-cube of lUnwii] {1983) ^nd recom- 
puted the rotation velocities by fitting a tilted-ring model to 
the velocity map. The rotcur task (Begeman 1989) of the 
GIPSY package (van der Hulst et al. 1992) was used for that 
purpose. No velocities are given inside a radius of 25', where 
the deficiency of neutral gas introduces large uncertainties, 
and where it is probable that there are large non-circular mo- 
tions. The error-bars for the inner RC are given by the differ- 
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ence between the velocities of the approaching and receding 
sides of the galaxy, except for the very few points where the 
formal error given by rotcur is chosen because it exceeds this 
difference. 

The inner part of the new RC compares well with the com- 
posite one presented in Widrow et al. (2003), which was com- 
piled from several data sets. The curves of the receding and 
approaching halves are shown to illustrate the symmetry of 
the gas motions in the galaxy: except, perhaps, near R ^ 70' 
where a difference of ^ 29 km s"' is observed, similar values 
are derived for the two sides of the galaxy. This asymmetry 
was also observed in previous Hi data (Braun 1991). 

The points ai R > 90' were computed from the new 
single-dish spectra, using the same inclination as for 
the inner parts (77°) a nd a systemic velocity of -300 
kms"' (de Vauc ouleurs e t al. 1991). The associated errors are 
taken to be equal to the largest error found in the inner part of 
the curve. The most striking feat ure of our ne w curve is that, 
far from having a declining RC JBraunlll991h . the velocities 
between 20 and 35 kpc are found to remain nearly constant at 
^ 226 kms"' . 

Models of the M3 1 warp at large radii imply a slightly in- 
creasing inclination and decreasing position angle of t he ma- 
jor axis as a function of radius ( Newton & Emersonj [T977t 
lBriggsllT990 ). As a consequence, use of a constant inclina- 
tion of 77° instead of a higher one at large R could lead to 
an overestimate of the rotation velocities, but only at most 
2.5% (or 6 kms"' ) with respect to the purely edge-on case. 
As for the position angle, if the Hi spectra are not taken ex- 
actly along the major axis, the derived velocities would be an 
underestimate of the true rotational velocity. This strength- 
ens our conclusion that the RC is not declining. Finally, the 
adopted error of ±29 km s~' is very conservative, and likely 
an overestimate of the errors which might arise, e.g., from a 
bad choice of the orientation parameters. 

Braun's conclusion that the M31 RC was declining at 7? > 
100' was probably determined by his two points at the greatest 
R, which were measured very close to the galaxy's minor axis 
(Tab. 3 and Fig. 8 of Braun 1991). For such an highly inclined 
galaxy, the outer regions close to the minor axis suffer from 
large deprojection errors. Furthermore, the superposition of 
several gas orbits along the line-of-sight due to the external 
warp is very important close the minor axis. These effects 
may lead to a bad estimate of the rotational velocity for the 
last two points of Braun's RC. A constant velocity at large 
radii, as the trend we derive, is in agreemen t with the flat RC 
model proposed bv lB rinks & BurtonI Jl984h . 

4. MASS DISTRIBUTION OF M3 1 

In order to model the mass distribution of M3 1 using the 
RC shown in Fig. |2l the luminosity profiles of the luminous 
components (bulge & disk) were derived from the B-band pa- 
rameters given in Walterbos & Kennicutt (1987, 1988). The 
Hi radial profile of Sofue & Kato (1981) was used for the 
gaseous disk, multiplied by a factor of 4/3 to account for He- 
lium. The observed colors (0.8 <(B-V)< 1.0), combined with 
population models, lead t o a value for the stellar component 
of 2.8 < (M/Lb) < 6.5 ( B ell &de Jong 2001 ). 

Figure |2l shows the best-fit mass model for M3 1 . The 
model is described in Carignan (1985), and has a dark mat- 
ter (DM) component represented by an isothermal sphere. A 
(M/Lb) of 7.5 (in units of M0/L0) is derived for the stel- 
lar disk. This is somewhat higher than the value expected 
from the colors, which means that this can be considered as 
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Fig. 2. — Rotation curve and mass model for M31. The new rotation 
velocities from the Effelsberg and GBT 100-m observations are for S > 21 
kpc. The velocities for R < 21 kpc are recomputed from the Unwin 1 198^) 
Hi data. The upper pointing triangles (light grey) are for the receding side 
and the bottom pointing (dark grey) ones for the approaching side, as obtained 
from a tilted-ring model (see text). The solid line is the best fit to the data. 
Each mass component is identified. A mass of 1.0 X 10^ Mq is used for the 
central black hole (Bender et al. 2005). 

a lower limit for the DM component (maximum disk model). 
Such large (M/Lb) values from bes t-fit mod els are not un- 
common: for example , value s of 9.4 ( lBegeman.1989.) and 8.5 
(|Blais-Ouellette et al.1 l 1999ft have been found for the galaxy 
NGC 3 198. Comparable estimates of (M/Lg) have been found 
from measurements of the vertical stellar velocity dispersion 
in the nearly face-on galaxies NGC 628 and NGC 1566 (van 
der Kruit & Freeman 1984). 

For 7? < 35 kpc, the model gives a stellar mass Msteiiar = 
2.3 X 10" Mq , an Hi mass Mh!= 5.0 x 10"^ Mq , and a 
dark matter mass of A4dm = 1.1 x 10" Mq , leading to a 
total mass of Mtot = 3.4 x 10" Mq . This translates in a 
MdarklMium ^ 0.5 and a total dynamical mass-to-light ratio 
(M/L);/vn ^ 6 at the last observed velocity point. 

5. CONCLUSIONS AND FUTURE WORK 

The main conclusions of this Letter are: 

• The present study has extended the measured rotation 
curve of M31 out to ^ 35 kpc. Contrary to previous 
studies (Braun 1991), the RC does not decline steadily 
from the center out to the last velocity point but remains 
nearly constant at ^ 226 km s"' for 20 < < 35 kpc. 

• A total mass for M31 of 3.4 x 10" MqIS derived 
for R <T>5 kpc. This is very similar to the mass of 

2.8 X 10" Al0(for R < 3\ kpc) found using kine- 
matical data of planetary nebulae ( Evans & Wilkinso^ 
EpOO), and can be compared with the MW mass of 

4.9 ±1.1 X 10" Mq within 50 kpc (Kochanek 1996) . 

• A study of the mass distribution gives a lower limit to 
the ratio MdarklMium ^ 0.5 at the last observed veloc- 
ity point. This is quite different from the no Dark Mat- 
ter model of Braun (1991). 

This study is the first part of a program to derive the best 
possible rotation curve for M31. Our observations confirm 
that Hi can still be detected past 30 kpc and that the outer 
RC is flat and not declining. An ideal RC will combine the 
following two data sets: 

(1) High resolution optical 3-D Ha kinematical data. 
These data will be obtained with the Fabry-Perot system 
FaNTOmM ( Gach et al...2002: .Hernandez et al. 2003„) . With 
the ^ 14 X 14 arcmin^ field-of-view at the Mont Megantic 
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Observatory telescope, it will be possible to get a mosaic of 
the whole optical disk with ~ 20 fields. This will allow us to 
get sufficient spatial resolution in the crucial rising part of the 
RC (|Blais-Ouellette et al. 1999) and at the same time provide 
the 2-D coverage necessary to derive properly the kinematical 
parameters. 

(2) Arc-minute resolution, high sensitivity synthesis 
Hi data. Because of the lack of short spacings, the VLA ob- 
servations (Braun 1990) are not sensitive to the large scale 
Hi structures expected in M3 1 . However, the DRAO synthe- 
sis array is sensitive to such structures and retrieves most of 
the Hi present in M3 1 . This has proved to be t he case for very 
extended nearby systems such as NGC 6946 dCari^nan et al.l 



[T990I) and DDO 154 jCarignan & Beaulieull98^ . 

These observations are in progress and should provide the 
optimal kinematical information for the best possible study of 
the mass distribution in M3 1 . 
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